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Abstract. This paper is devoted to rigidity results for some elliptic PDEs and related interpolation 
inequalities of Sobolev type on smooth compact connected Ricmannian manifolds without boundaries. 
£f) ' Rigidity means that the PDE has no other solution than the constant one at least when a parameter is 

in a certain range. This parameter can be used as an estimate for the best constant in the corresponding 
f"*"^ , interpolation inequality. Our approach relies in a nonlinear flow of porous medium / fast diffusion type 

■ which gives a clear-cut interpretation of technical choices of exponents done in earlier works. We also 

' establish two integral criteria for rigidity that improve upon known, pointwise conditions, and hold for 

' general manifolds without positivity conditions on the curvature. Using the flow, we are also able to 

rZr , discuss the optimality of the corresponding constant in the interpolation inequalities. 
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In the past decades there has been considerable activity in establishing sharp inequalities using maps or 
flows. The basic idea is to look for a flow on a function space along which a given functional converges to its 
optimal value, i.e., one turns the idea of a Lyapunov function, known from dynamical systems theory, on 
its head. An example is furnished by the relatively recent proofs of the Brascamp-Licb inequalities using 
non- linear heat flows in [171 1571 126] . Using the same methods a new Brascamp-Lieb - type inequality on 
S d was proved in [55] . Likewise the reverse Brascamp-Lieb inequalities can also be obtained in this fashion 
| (see [2]). Another example is the proof of Lieb's sharp Hardy-Littlewood-Sobolev inequality given in [2"7] 

where a discrete map on a function space was constructed whose iterations drives the Hardy-Littlewood- 
Sobolev functional to its sharp value. Likewise, the sharp form of the Gagliardo-Nirenberg inequalities due 
to Dolbeault and Del Pino can be derived using the porous media flow (see [29] and [25] for the relation 
between the porous media equation and a special class of Hardy-Littlewood-Sobolev inequalities). 

! Closely related are the proofs of sharp inequalities using transportation theory. The earliest use of 

£f) • transportation theory to our knowledge was in Barthe's proof of the Brascamp-Lieb inequalities as well 

as their converse, in |13j . Transportation ideas were also applied in |32j for proving the sharp Gagliardo- 
Nirenberg inequalities and in [S3] for proving sharp trace inequalities for Sobolev functions. 

In this paper we invent a new type of porous media flow on Ricmannian manifolds that allow us to give 
relatively straightforward proofs as well as generalizations of rigidity results of [3T] [TTJ 25] for a class of 
d • non-linear equations. Before describing the flow, wc discuss the rigidity results. 

Throughout the paper we assume that (DJl, g) is a smooth compact connected Riemannian manifold of 
dimension d > 1, without boundary with A g being the Laplacc-Beltrami operator on 971. For simplicity, 
we assume that the volume of 9Jt, vol(9Jt) is 1 and we denote by dv g the volume element. We shall also 
denote by D\ the Ricci tensor. Let Ai be the lowest positive eigenvalue of — A g . We shall use the notation 
2* := ^ if d > 3, and 2* := oo if d = 1 or 2. 

Let us start with results dealing with manifold with curvature bounded from below and define 

p:=inf inf 9t(£,f) 

Theorem 1. Let d > 2 be an integer and assume that p is positive. If X is a positive parameter such that 
A < (1 — 0) Ai + — whe~~ 



X 
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d-1 d(d + 2)+p-l 
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then for any p € (2, 2*) ; £/ie equation 
(1) 



p-2 



has a unique positive solution v 6 C 2 (9Jt), which is constant and equal to 1. 

Such a rigidity result has been established in [TTJ Inequality (1.11)] by D. Bakry and M. Ledoux (also 
see [39]), and in (49] Theorem 2.1] by J.R. Licois and L. Veron. An earlier version of it can be found in a 
work of M.-F. Bidaut- Veron and L. Veron, |21j . with an estimate only in terms of p. which is not as good 
for general manifolds but coincides with the one given above for spheres. Earlier results can be found in 
|40j . Each of these contributions relies either on the Bochncr-Lichncrovicz-Weitzcnbock formula or on the 
carre du champ method, which are equivalent in the present case. 

The case of the critical exponent p = 2* when d > 3 is not covered by Theorem [JJ and requires more 
care, due to compactness issues and conformal invariance. For simplicity we shall restrict ourselves to the 
subcritical case in this paper, at least as far as rigidity results are concerned. 

We shall now state a slightly more general result, which does not require that p is positive but only 
involves the constant 

(1 - 6>) (A g u) 2 + -p^r £R(Vu, Vu)l dv g 



inf 

mGH 2 (OT) 



an 



d- 1 



Here 9 is defined as in Theorem [TJ If p is positive, it is not difficult to check (see Lemma [7]) that 

A* > (l-9)X 1+ 9-^- . 

a — 1 

Here we shall simply assume that A* is positive but make no assumption on the sign (in the sense of 
quadratic forms) of the Ricci tensor 9t. Our first result generalizes Theorem [JJ as follows. 

Theorem 2. With the above notations, if X is such that 

< A < A* , 

then for any p G (1, 2) U (2, 2*), Equation (TTJ) has a unique positive solution in C' 2 (VJl), which is constant 
and equal to 1. 

In this statement we include the case p € (1,2). We may observe that lim p ^i + 6(p) = 0, so that 
limp-*.!, \*{p) — Ai, as soon as p is bounded (but eventually negative). Notice that if SEJl = E> d , then 
A* = Ai = dp/{d — 1) = d, p = d — 1. In this case, the result of Theorem[T]is then optimal, but this was 
already known from [211 Theorem 6.1]. As we shall see in Proposition rj31 we always have 

A* < Ai 

and (TTJ) has non constant solutions for any A > Ai. We can actually give another integral criterion for 
rigidity that slightly improves the result of Theorem [2] Recalling that 9 defined in Theorem Q] is given by 

(d-lf{p-l) 



9 = 



let us define 



UgU 



9 A 

-AgU 



where H g u denotes Hessian of u, and define 



(2) 



A* := inf 

«GH 2 (Olt)\{0} 



d(d + 2)+p- 

(d-l)(p-l) 
9 (d + 3 - p) 

d 



Vu <E> Vu 



(l-fl) / {AgUYdVg + 



an 



an 



9 IV^I 2 
d u 



$H(Vm, Vu) 



\Vu\ 2 di 



an 



It is obvious that 

A* < A, , 

and it is not hard to see that there is equality in case of the sphere. 
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Theorem 3. For any p S (1,2) U (2,2*), Equation {T]) has a unique positive solution in C 2 (2t) if A € 
(0,A*), which is constant and equal to 1. 

Our approach relies on a nonlinear diffusion equation of porous media (fast diffusion) type. Its interest 
is that it simplifies the proofs of previously known results. It applies to equations which are not of power 
law type (see Theorem I14[) and also provides an expression of the error term (see Corollary I15[) in the 
interpolation inequality 

(3) \\Vv\\l* (m > -A- \\\v\\l P(m) - \\v\\ 2 h2m) ] V.Gtffl). 



p-2 



Notice that this inequality makes sense not only for p <E (2, 2*) but also for p € (1, 2) because in this last 
case we have |H|Lp(9Jt) < MIl 2 (ok) by Holder's inequality and p — 2 < 0. Inequality © has been stated 
in [3TJ Corollary 6.2] with A = dp/(d — 1), but has also been established in the case of the sphere in [TB] 
with A = d using spectral estimates and Lieb's sharp Hardy-Littlewood-Sobolev inequality [50]. We shall 
see in Proposition [13] that A* < Ai. As a consequence of Theorem [3] we get the following result. 

Theorem 4. For any p e [1,2) U (2,2*] if d > 3, p G (1,2) U (2, oo) if d = 1 or 2, Inequality © holds 
with A = A 6 [A*, Ai]. Moreover, if A* < Ai, then the best constant A is such that 

A* < A < Ai . 

If d = 1, then A = Ai. 

While the case p = 2* is not covered in Theorem [3] inequality ([3]) holds in that case by taking the limit 
p — > 2*_ when d > 3. In the limit case p = 1, inequality ([3]) is nothing else but the Poincare inequality 
which is consistent with the fact that lim p _>i + A*(p) = lim p _^i + A*(p) = Ai. Using u = l+ £(/?asa test 
function where ip is an eigcnfunction associated with the lowest positive eigenvalue of — A g , it is a classical 
result that the optimal constant in ([3J is A < Ai. For the same reason, a minimum of 

A 



v i ^ \\Vv\\l 2m) 



P 



MlL(OT) - IMlL=(0Jt) 



under the constraint |M|lp(ot) = 1 is negative if A > Ai . Incidentally, such a minimization provides a non 
constant nonnegative solution of ([1} (see (23] [37]). It is natural to wonder what happens when p = 2. The 
reader is invited to check that by taking the limit as p — > 2, one gets a logarithmic Sobolev inequality, 
which is of independent interest (see for instance [21] or [37] for more results in this direction and further 
references) . 

The main contribution of this paper is to introduce the flow 

(4) Ut=u *-*pU gU + K \^ 

where we set 

(5) k = 1 + P (p - 2) . 

The function u which evolves according to ([4]) and the function v which appears in Theorems [2] [3J and 
Theorem [4] are related by v = u 13 . The flow (J4| contracts both sides of Inequality ([3]) and also 

2/p- 

(6) T[u]:= I \V{u p )\ 2 dv g + ^: / u 2p dv n -[ / u Pp dv„ 



A 

p-2 



9 

wi Van 



This contraction property actually determines j3 in a way which is somewhat more transparent than in 
[21] [49] (see Section[3]) and allows for simple generalizations, like the one stated in Theorem IT4l concerning 
more general nonlinearities than just power laws. The choice of k is determined by the condition that 
|M|lp(ot) ^ s invariant under the action of the flow (see Lemma [SJ). The equation ([4]) is of the type of a 
porous media equation. Let us note in this context that Hamilton's Yamabe flow can be seen as the fast 
diffusion case of a porous media equation |59] but it is not clear whether there is any connection between 
this and 

This paper is organized as follows. In Section [2] we introduce two simple computations, which differ 
from the ones that can be found in (21] [11] [49] : results are stated in Lemmata [5] and [6] Section [3] is devoted 
to the evolution of the functional T given by ([6]) under the action of the flow defined by ([4]). Proofs of 
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Theorem [2] and Theorem [4] arc given in Section [4] However, a simple inspection of the proofs allows to 
state an improved version of Theorem 2] with an integral remainder term. The corresponding results and 
some consequences have been collected in Section[5j as well as the proof of Theorem[3] The question of the 
optimality will also be discussed in Section [5] The use of the flow ((4]) shows that estimates are in general 
not optimal, although our approach relies only on global estimates (rather than on pointwise estimates). 
We stress the fact that global estimates allow us to consider the case of sign changing curvatures. 

Before dealing with our method, let us give a brief list of papers which are relevant for our approach. 
We will only mention the ones which are closely related to our purpose and select in the huge literature 
corresponding to the various topics listed below the earliest contributions we are aware of, or some recent 
papers to which we refer for more details on the historical developments. The motivation for this work 
comes from the method developed by B. Gidas and J. Spruck in [40] Theorems B.l and B.2] for proving 
rigidity and later improved first by M.-F. Bidaut- Veron and L. Veron in [2l] Theorem 6.1], and then by 
D. Bakry and M. Ledoux in [XTJ Inequality (1.11)] and by J.R. Licois and L. Veron in [49] Theorem 2.1]. 
The choice of the exponents was somewhat mysterious in these papers while the choice of (3 and k is rather 
straightforward and natural in our setting. 

When 9Jt = S d , the rigidity result with optimal range in terms of the parameter A was established in [21] 
and later, by other methods, in [16] ■ The case of the critical exponent was known much earlier since the 
inequality corresponds to the Sobolcv inequality on the Euclidean space, after a stercographic projection: 
we can refer to [4] Corollaire 4], to |54) for issues dealing with conformal invariance, to |36j for results on 
the sphere, and to |37j for consequences in the spectral theory of Schrodingcr operators. 

For general manifolds, the so-called A-B problem, which amounts to determine the two best constants 
in interpolation inequalities, has attracted the attention of many authors and it is out of the scope of 
this paper to list all of them. Let us simply quote [5j [6] [47] [44] [43j [12] as a selection of contributions 
on the "second best constant" problem. The interested reader can also refer to [42] for a large overview 
of the topic. With some exceptions (see for instance f24j [35]), almost all related papers are concerned 
either with the case of critical exponents and related concentration issues or with the limiting case p = 2, 
where Inequality ^ has to be replaced by the logarithmic Sobolev inequality (see [35] [UJ [5S] E2] for early 
contributions to this topic). 

Inequality (j3j has been studied for p e (1,2) and various probability measures, mostly defined on the 
Euclidean space: we may refer for instance to [HI QUI EES EM (SB [22] ■ The paper [35] by P. Deng and 
F. Wang is of particular interest, as it gives estimates of the constant A in ([3]) for arbitrary p 6 (1,2) 
using several methods, including the case of compact manifolds without assuming the positivity of the 
curvature. In the particular case of the sphere, the whole range p G (1, 2*] was covered in the case of the 
ultraspherical operator in [HI [19] . Also see [36] for the use of the ultraspherical operator in the context 
of rigidity results. 

Using flows on manifolds is a classical idea, including flows associated either to the heat equation or 
to porous media equations, or in relation with Ricci (see for instance |31j ) and Yamabe flows (see |58| 
Appendix III] for an introduction). Another line of thought based on flows refers to the carre du champ 
method: see [HI El HI IB] an d 0S! for a quite detailed review by M. Ledoux. In the framework of diffusion 
equations from the PDE point of view we primarily refer to [3] in the linear case, and to [29] [33] [28] for fast 
diffusion equations in the Euclidean space. This approach was later interpreted in terms of gradient flows 
with respect to Wasserstcin's distance in |55j , and then extended or reinterpreted in various contributions 
in the context of entropy methods. Such a direction of research is out of the scope of the present paper, 
although it has been a strong source of inspiration. Much less is known in the case of Riemannian 
manifolds: see |341 [5Tj . As far as we know, the flow defined by (f4]) has not been considered before, except 
in [36| when /? = 1, or when k = and either p = 2 of p = 2*. 



2. Some preliminary computations 
If Ajj and B-jj are two tensors we use the notation 

A • B := g l > m g>' n A ld B m> „ and || A|| 2 := A • A . 
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Here g 1,3 is the inverse of the metric tensor, i.e., g 1,3 gj.k = S l k where we used the Einstein summation 
convention and 8\. denotes the Kronecker symbol. Denote by H s u the Hessian of u and by 



L g u := H g u — — A g u 



the trace free Hessian. The following lemma is well known, but since it is of some importance for what 
follows we give a short proof. 



Lemma 5. If d > 2 and u 6 C 2 (SOT), then we have 



(A g u) 2 dv g 



an 



d- 1 



an 



L g u\\ 2 dvg + j— j 



9\{Vu,Vu)dv. 



an 



9 ' 



Pivof. Integrating the Bochner-Lichnerovicz-Weitzcnbock formula 

1 
2 

yields 



an 



A |Vu| 2 = ||H s u|| 2 + V(A g u) • Vu + JH(Vu, Vu) 
(A 9 u) 2 dw 3 = / ||H g u|| 2 dv 9 + / £H(Vu, Vu) 



an 



an 



where, incidentally, we may remark that the right-hand side is nothing else than the integral of T2(u,u) 
in the carre du champ method: see for instance Next we notice that 



[L g u] -3 = 0, 



because A g u is the trace of H ff u, which implies 



an 



hgU \\ 2 dVg 



an 



\n g u\\ 2 dvg - - 



an 



{AgUfdVg , 



from which Lemma [5] follows. 



□ 



The next lemma will be useful to express some important quantities as a square. Here we deal with the 
cross term. 



Lemma 6. If u € C 2 (9Jt) is a positive function, then we have 



an 



a |V^| 2 , d 
A g u— r dv g = J — 



an 



|Vu| 4 , 2d 
— dV °-—2 



ir 



[ Ml- 


'Vu (8 Vu" 




/an 


u 



rfu„ 



Proof. An integration by parts shows that 

|Vu 12 



an 



A ff u- 



■ dv n 



|Vu| 4 , 
an u 



/ [H s u]- 


'Vu ® Vu" 




/art 


u 



dv 



9 ■■ 



which equals 



an 



|Vu| 4 , 
u 2 



Jan 



Vu ® Vu 



an 



a |V«| 2 

A u dv a 

u 



by definition of L„u. Lemma [5] follows. 



□ 



Recall that g) is a smooth compact connected Riemannian manifold without boundary and consider 
the Poincare inequality 



(7) 



/ |Vu| 2 dv g >Ai / |u-u| 2 du g VueH^SDt) such that u= udv g , 
Jm Jm Jm 



with optimal constant Ai > 0. By standard compactness results, it is known that Ai can be characterized 
as the minimum of |Vu| 2 dv g j J w \u — u| 2 dv g on all nonconstant functions u £ H 1 (9Jl) and is achieved 
by a function <p S H 1 (93T) such that ||</5||L 2 (art) = 1 an d ~<p — 0. 

The third lemma of this section is a standard result, that is nonetheless crucial in our approach. For 
completeness, let us give a statement and a short proof. 
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Lemma 7. With the above notations, we have 



[ (A g u) 2 dvg > Ai / \\7u\ 2 dv g V«eH 2 (9J!) 

JWl J 931 



Im J 931 

Moreover, Ai is t/ie optimal constant in the above inequality. 

Proof. A simple computation based on the Cauchy-Schwarz inequality 



2 / r \ 2 

|2 j„, \ _ I A „. tt\ J„. \ ^ I f \ „.\2 j„, / L. — 12 . 



shows that 



/ |Vu| 2 <ivg) =1 / A g u (u — u) di> g I < / (A g u) z dv g I \u — u\ z dv g 
Jm ) \Jm ) Jm ' Jm 



(AgUfdVg> ^ U? *f 

Jm \u-u\ z dvg 



and one observes using ([7]) that 



(f m \X7u\*dv,:'- 



Jot \u-u\ 2 dv 



> Ai / |Vu| 2 dv 9 , 



m | U - u| i« s J 931 

which concludes the proof. Notice that optimality in the Cauchy-Schwarz inequality means that A s m and 
(u — u) are collinear, that is, (u — u) is an eigenfunction. It is then clear that equality holds if u = ip, an 
eigenfunction associated with the first positive eigenvalue Ai. □ 

As a consequence of Lemmata [5] and [3 we find that 

Ai = inf Im ( A 9 U ) 2 dv 9 > d ^ Jgg g(Vu, Vu) rf^ g 



J OT |Vu| 2 dwg d-l J m \\7u\ 2 dvg 

thus establishing the Lichnerowicz estimate: Ai > -pj (see [3H1 P- 135] for the original statement, and pZOl 
p. 179] for the Obata-Lichnerowicz theorem). 

3. A POROUS MEDIA FLOW 

This section is devoted to the study of the functional defined in © and its evolution under the action of 
the flow defined by ([3]), which is an equation of porous media type (the exponent /3 is in some cases smaller 
than 1 and one should then speak of fast diffusion: see discussion at the end of this section). We will 
expose the computations without taking care of regularity issues that are requested to justify integrations 
by parts and give afterwards a sketch of the proofs and regularizations that are eventually needed. The 
first statement explains how the coefficient k is chosen. 

Lemma 8. With n given by §5§, the functional u t— > u@ p dv g remains constant if u is a smooth positive 
solution of 

Proof. Computing the time derivative, we get that 



dt 



an 



Jm V 



dvg = 



by the choice of n. □ 

Next we compute the time derivative of the remaining terms in the functional 
Lemma 9. If if u is a smooth positive solution of (j4|), the following identity holds 



1 d t\ ^ 



931 



(AgU) 2 + ( K + P-l)AgU l -^ + K(P-l) 



u 



dvg + A / | Vu| 2 dv g . 
m 
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Proof. By Lemma [8l we only have to compute ^ J OT |V(w' 3 )| 2 dv g and J m u 2 ^ dv g . We get 



2 A 



|V(u")| a dv fl 



an 



|Vu| : 



an 



£ 2 / A fl «+G8-1) 



|Vm| : 



|Vu| : 



/3 2 



an 



|Vu| : 



(Agit)- 4 + (« + p - 1) A s w J — + k(/3-1) 



and 



1 d_ 

2 dt 



\Vu\ 



dv„ 



(8) 



= /? / u 2 / J - 1 « a - 2 ^ (a 3 u + 
= P(k-1) \Vu\ 2 dv g 
= P 2 (p-2) f \Vu\ 2 dv g , 

which completes the proof. 

Notice that for any 9 S (0, 1) we can write the result of Lemma [9] as 

^ j/W = (1 - 0) J (A 3 «) 2 dv g G[u} + A £ \Vu\ 2 dv 



□ 



where 



a.n 



V?;l 2 

(A 9 u) 2 + (« + /3 - 1) A s u J L + „ (p - i) 

U M z 



|Vu| 



9 ' 



We shall now introduce a key quantity for our method. For any positive function u <G C 2 (9Jt) set 

~Vu®Vu g |Vu| 2 " 



Q«u:=L fl u-J^-|(K + )8-l) 

Notice that is a matrix valued, traceless quantity. The exponent /? in the flow (Q| and the value of 
9 G [0, 1] will be chosen below in this section. 

Lemma 10. Assume that d > 2. With the above notations, any positive function u € C 2 (5D?) satisfies the 
identity 

rl 

m{Vu, vu) dv g 

'an 



G[u] = 



d-1 



lull 2 dv a 



an 



an 



|Vu| 4 , 

T~ dv 9 

u 



W ith fx := I + /3 - l) 2 - k (/? - 1) - (« + P - 1) 
Proof. By Lemma [5] and [6l we know that 



0d 
d-1 



an 



L ff u || 2 dw s 



<K(Vu, Vu)di 



+ (/s + 0-l) 



_d + 2 

which can be arranged into 
9d 



'a.n 

|Vu| 4 



a.n 



2d 



/ Ml- 






+ K (P-1) [ 


dw 9 




u 







|Vu| 4 
— 5— dv s 
u 2 



0M = 



d-i 



an 



T 112 j 2d-l 



[ M- 


'Vu g) Vu" 


dv 3 


/on 


u 





+ «(£ -!) + (« + P-1) 



d + 2 



an 



|Vu| 4 9d 

5~~ dv 9 + 3 7 

u 2 d-1 



9t(Vti, Vu) du 9 



an 



Noting again that [L s u] -5 = 0, we can rewrite the first line as 



9d 
d-1 



L„u II dii„ 



2 d- 1 



a.n 



d + 2 



K + 0-1) [ [LgU] 



Vu<g> Vu g |Vu| 



d u 



dv„ 
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which, by completing the square, turns into 



Od 
d-1 



an 



WQtufdVg- 



1 



1 



e \d + 2 



|Vw| 4 , 
on u 



because 



an 



Vu® Vu g |Vu| 2 
u d u 



dv„ 



|V U | 4 
an u 



This completes the proof of Lemma [TQl 



□ 



Now let us explain how j3 and 9 arc chosen. Up to now, we have shown that 
'an a — j- ijfjji - jcqi 



\\Q e q u\\ 2 dv g 



5H(Vu, Vu) dt; s 



/' 



|Vu| 



an 



di; s + A / |Vu| dv fl 



an 



where, after taking into account ([5]), /i is given by 



d + 2 



d+2 



(p-l) 2 -(p-2) /? 2 



d + 2 H 



1 



The goal is to find values of /3 for which is equal to zero, thus maximizing the coefficients of the other 
terms. Unless 6 (p - 2) = (p - l) 2 (d - l) 2 /(d + 2) 2 , the coefficient /i is quadratic in terms of /3, and we 
can choose the extremum of /i with respect to /3, i.e., 

(d + 2)(d + 3-p)6> 
(d - l) 2 (p - l) 2 - (d + 2) 2 (p - 2) ' 

With this specific value of (3, we adjust the value of 9 so that /z(/3) = 0. All computations done, we find 
that 

w d(d + 2)+p-l p d+3-p 

and one can easily check that with these values, and without condition on 9, we actually have /.t = 0. 
When p varies in the interval (0,2*), 9 = 9(p) ranges from 9(1) = to lim p ^2* 0(p) = 1 if d > 2. It is 
left to the reader to check that the expression of Q g given in Section [1] coincides with for the above 
choices of 9 and f$. Notice that with such a choice of 9, we get (3 > 1 for any p G (1, 2*) if d > 3, and any 
p G (1, 5) if d = 2. On the contrary, we find that j3 < if p > 5 and d = 2. If p = 5 and d = 2, then we 
have p = (i — 3) /3 2 + 1, so that 9 has to be chosen strictly less than 1/3. 

Collecting the results of Section [3j with A* given by ([2]), we have found that 

Proposition 11. Let d > 2, p 6 (1,2) U (2,2*) and assume that (3 and 9 are given by (|9|) if p ^ 5 or 
If p = 5 and d = 2, we assume that 9 e (1/3, 1) and /3 = ^9/(39 — 1). MiA re given 6y ((S|), we 

#ei 

2^1^<(A-A /jV.| 2 d, g 
if u is a smooth positive solution of Q • 

The case d = 1 has to be handled separately because Lemma [5] does not make sense. However, 
computations are much easier because 



dv„ = 



1 



an 



,'|4 



3 7 OT u 2 



dw„ 



so that there is no cross term to compute. The reader is invited to check that Lemmata [JJ [5] and [5] hold 
when d = 1, thus showing that 



1 ^TTf 1 

2^ dt™ 



an 



dv g + A / \u | di> 9 . 
an 
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with n := \ (k + fi - 1) + K{fi - 1) = 3 P {p - 1) + (l + p (p - 2)) (P - 1). Hence, for all p > 1, there are 
two values of /3, given by 



(10) P ± := PZ1±V^EW±1 

V ' H± 3(p-2) 

such that /i = 0. For P = p±, we have found that 



= -J W'\ 2 + a W\ 2 dv g < 



for any A < Ai- 



Proposition 12. Let d = 1, p € (1, 2) U (2, +oo). Assume that k and p are given respectively by ([5]) and 
by ITU]). T/ien 



^^«]<(A-Ao/jV»|»*,. 
if u is a smooth positive solution of Q and 

(i) either p > 2 and /3 g (— cxd, p_] U [/?+, +oo), 

(ii) orpE (1,2) and /3 G [/?_, 0) U (0, ^+]. 

In order to justify the above computations, we have to establish that the solution of ^ is smooth and 
converges to a constant. Since this is not at all the scope of this paper, let us simply sketch that strategy 
of proof. A positive solution of ([U of class C 2 is actually smooth by standard elliptic theory, and thus 
we know that the initial datum for (|4]) is smooth. If we simply consider an optimal function for (|3]), wc 
may proceed by regularization and approximate the solution by a bounded solution, which is also positive, 
bounded away from zero. A further regularization of the initial datum allows to get back to the setting of 
Theorems Q][3l By applying the Maximum Principle, we can then extend the boundedncss properties (from 
above and from below) of the initial datum to the solution of (UJ) for any positive time. Regularity results 
of Harnack and Aronson-Bcnilan type like the ones in |51j and subsequent papers then hold. Higher order 
regularity results follow: see for instance [58] in the Euclidean case. The delayed regularity phenomenon 
can be avoided by starting with smooth enough initial data. 

It remains to show that the solution converges to a constant. One can argue as follows. Using Lemma[8] 
and Holder's inequality, we know that J m u 2 @ dv g is bounded uniformly in t from above if p > 2 and from 
below if p G (1,2). By (|8|), is bounded uniformly in t and, as a consequence, J + °° |Vn| 2 du 9 ) dt is 
finite. Standard arguments allow then to prove that lim t ^ +oc j m |Vu| 2 dv g = 0. 

4. Proofs 

Let us consider A* and 9 as defined by @ and 
Proposition 13. Using the notations of Section^ we have the estimates 

(1-0)\ 1+ 0J*£- <A*<A*<Ai . 
d — 1 



Proof. The first inequality is a consequence of Lemma [7J For the second one, we observe that 

1 

e^O e 2 

Using ti=l + £»asa test function, we find that 



Km^\\Ql(l + 6v)\\ 2 = \\L g v\\ 2 . 



A* < inf Y^L^ 3 - *i 
«eH>(0K)\{0} f m \Vv\ 2 dv g 

according to Lemma [5] and Lemma [7] □ 



Proof of Theorem If v — u@ is a solution to (fTJ) , it is straightforward to check that 

|Vu| 2 A 

0(p-2) 



A 9U +(/3-l)J L_ ^ (u-u K )=Q 
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Using the fact that J m u K (^A g u + k J dv g = 0, it follows by Proposition [TT] that at t = wc have 



= - 



an 



A g u + 08-1) 



A 



|Vd : 



u p (p - 2) 



|Vu| : 



dv a 



so that u is constant if A < A*. 



□ 



Proof of Theorem^ We apply Proposition 1111 to u(t = 0, •) = v 1 ^. For A < A*, J-[u] is nonincreasing. 
We know that as t —> oo u and Vu converge respectively to a constant and to 0, which proves that 
J-[u{t = 0, •)] > J-[u(t, ■)] > = limt^oo ^[uit, •)] for any t > 0. Hence, for any A < A*, this proves that 
J-[u] is nonnegative for any u, i.e., that © holds. 

If A = A* and if v = is a solution to (UJ, we get that 



< T[u] = \\Vv\\l 2 



A. 



(an) 



p-2 



IMIlp(ot) - INIi^e 



(SB) 



A, 

p-2 



II W IIlp(9JI) _ IMIlp(<&1) 



thus proving that |H|lj>(ott) ^ 1- As a special case, we get that if v is an optimal function for ([3]), then we 
can assume that ||u||lp(s»i) = 1 with no restriction and if u is the solution of @ with initial datum u 1 /' 3 , 
then 



) I (A g u) 2 dv g + -^- 
an a - i 



an 



\q gU \\ 2 + m(Vu,Vu) 



dv Q — A* / [Vitl di>,; 



an 



. 



Arguing as in Proposition [T3l since ||Q S 



|L 9 u|| 2 as t — > +oo and using Lemma [SJ we find that 



0= (1-6) f [A g u) 2 dv g + -j—— 
Jm a — l 

~[(l-0) / (AgU) 2 dVg + 



an 



\Q g u\\ 2 + d\{Vu,Vu) 



dv q — A* / |Vit| du„ 



an 



an 



|L s w|| 2 + $H(Vm,Vu) 



dw — A* / I Vitl 2 c?w„ 



an 



an 



(A g u) 2 dvg - A* / |Vu| 2 du s > (Ax - A*) / |Vu| 2 d?; 



an 



an 



thus obtaining an obvious contradiction if A* < Ax, unless u is constant. Hence, in that case, we know 
that all minimizers of T are constants when A = A*. 

Now let us consider the optimal constant A in If A = A*, consider a sequence (A„)„ e N such that 
A„ > A* for any n € N, lim„_ ) . 00 A„ = A* and a sequence (w n )nGN of functions in H 1 (971) such that 
ll t 'n||L 2 (an) = 1 f° r an Y weN and J-\ n [v n ] < where, with obvious notations, J-\ n denotes the functional 
T written for A = A n . Actually we can consider minimizers of T\ n in {u G H 1 (SDt) : ||u||l 2 (9ji) = !}• Let 
us consider the quotient 



Q[v] := 



(p-2)\\Vv\\l 



2 (an) 



IMIL>(ot) II u IIl 2 (ot) 



and observe that the optimal constant A is such that 

A = inf Q[v] = lim Q[v n ] . 

ueH^ajQVfO} n^oo 

Unless p = 2d/(d — 2), d > 3, up to the extraction of a subsequence, (w n ) n eN converges to a limit which 
solves the Euler-Lagrange equation and is therefore constant. Because of the constraint ||t | n ||L 2 (an) = L 
the limit is 1. Let us define w n := e n (v n — 1) with e n := \\v n — 1||l 2 (3H) — > 0. A standard computation 
then shows that 



Q[v n ] = Q[l 



|Vw„|| 2 2 ( OT ) 



K|| 2 



L 2 (OTt) 



from which we deduce that A* > Ax. Hence A* < Ai means that A > A*, which concludes the proof. □ 
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5. Improved results and consequences 

Theorem [I] is a consequence of Proposition 1131 The expression of A* in Theorem [3] is not as simple as 
the one of A* in Theorem [2] but provides an improved condition for rigidity. Up to now, we have dealt only 
with power law nonlinearities. It is not difficult to generalize our results to nonlinearities that compare 
with power laws, for instance f(v) = |u| p_1 v + H 9_1 v for some q > p. Our first extension of Theorem[3] 
goes as follows. 

Theorem 14. Let f be a Lipschitz increasing function such that /(0) = and 

J—(f'(v)-(p-l)iM]< VoOa.e. 
P - 2 V v J 

Then under the same conditions as in Theorem^ for any A £ (0, A*), the equation 



A g v 



■ v = f(v) 



has a unique positive solution, which is equal to a constant c such that X = (p — 2) /(c)/c. 



As a consequence, we also have the inequality 

A 



L 2 (OT) 



p-2 



pF(v) dvc 



a.n 



2/p 



Him 



l 2 (an) 



>0 V v e H 1 (9Jt) 



for some A with the same properties as in Theorem [4j and F[v] := f Q f(s) ds, so that F(v) = - \v\ p in 
the case of power laws. 

Proof of Theorem \14\ The proof is the same as the one of Theorem [3l except that 



F(v) dv g 



F(u p )di 



is not preserved by the flow defined by ((4]). However, the functional 



F[u] := ||V(</)||£ 2(OT) -^ 



U m pF{vP)dv g f p 



is still non increasing because, with v = vfi , we have 



F(v) dv g 



p-2 dt J m ~ "~ u p-2 Jm L 
The proof then goes as in the power law case. 



(/3 + K -l)M -/?/» 



\\7u\ 2 dv g > . 



□ 



Our method actually provides an integral remainder term that can be computed using the flow, when 
A = A*. If u a smooth solution of (U]) with initial datum v 1 ^ , let 



K[v] := (3 2 



(1-d) I {A g u) 2 dv g + 



3Ji 



Od 



|Qg«|r + tH(Vu,Vu) dvg-k* \Vu\ z dv g 



A, 
P-2 7o 



pF(u f3 )dv i 



a.n 



an 



(/? + K -i)^-/?/V) 



|Vd 2 du„ 



Corollary 15. FFii/i i/ie same notations as above and A* defined by for any p 6 (1, 2)U(2, 2*] i/d > 3 
and any p € (1, 2) U (2, oo) otherwise, for any smooth positive function v on 9Jt, we have the inequality 

, 2/p 

p F(v) dv g 



\Vv\\h(<m) 



A. 



p-2 



an 



MIl 2 (OT) 



> 2K[v] 



Let us conclude this section by some comments in the case p £ (1,2) and d > 2. By the definition © 
of A*, we know that 



A* > (1 - 9) Ai + 6 p+ with := 



inf 



f m \\Q g u\\ 2 + m(\7u,\7u) 



dv a 



d — 1 tieH 2 (an)\{o} 



J m \\7u\ 2 di 
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We also know that p* > p and hence A* > A*. Let us denote by A(p) the best constant in ([3]). From 
Theorem 2] and Proposition HU1 we have 

Ai > A(p) > A*(p) > (1 - 0) Ai + 

with 

(d-l) 2 r? 



d(d + 2) + ?? ' 
If p = 1, we have A(l) = Ai. 

On the other hand, passing to the limit as p — > 2 is straightforward so that we get 

|a(2) / M 2 log ( ^ ) < \\vv\\l* m v« e H^an) . 

z JOT V II v IIt ?«yrt\ / 



'an VII^IIl^ot) 



HereA(2) denotes the best constant and we have A(2) > A*(2) > (l-0)Ai+0p* with = ((d-l)/(d+l)) 2 . 
The interpolation method based on spectral estimates of [15] and later extended in [2] Theorem 2.4] (also 
see (36] Section 2.2]) shows that 

1-77 

Alp > Ai with a = X 1 /A(2) and r)=p-l . 

1 — ?/ Q 

Such an estimate is better than A(p) > (1 — 0) Ai + 6 p+ at least in a neighborhood of p = 2 if A(2) > 
(1 ^ 02) Ai + 02 p* with 02 = ((d — l)/(d + l)) 2 - This proves that for estimates on the best constant in ©, 
there is space for improvements. 

6. Concluding remarks 

All computations in this paper are based on the action of the flow ^ on the functional T . The 
functional T is monotonously decaying with respect to time for the right choice of the exponent /3 (and 
the optimal choice is therefore prescribed by the method) while the solution of (U]) converges to a constant. 
The method gives a meaningful way of reinterpreting the computations of [3TJ 3D] [HI EH] and suggest the 
improvements of Theorems [2] [3] and [4] However, the use of the flow itself is not mandatory except in the 
proof of Theorem [4] to handle the case A = A* and in the extensions corresponding to Theorem [14] and 
Corollary 1151 The flow approach however paves the road for the various improvements presented in this 
paper: the integral criteria defining A* and A* or the extension to general nonlinearities as in Theorem 1 141 
In Theorem 3] an alternative approach could be used to directly prove that 

/ (A g u) 2 dv g + [ [ ||Q ff u|| 2 + $H(Vu, Vu)l dvg - A* / \Vu\ 2 dv g =0 
Jot a— i J m L J J m 

means that u is a constant if A* < Ai. Such a result is a consequence of our method, at least when u is a 
solution of (U) whose initial datum is a minimizer of T . 

The choice ^ is the one that minimizes the value of and therefore maximizes the value of A*. 
However, if we relax a little bit the condition on £ (0, 1] and take such that 

(d-lf{p-l) 
d(d + 2)+p-l ' 

then there is a whole interval [/3_,/3+] of values of /3 for which the method of Section [3] still works, i.e., 
for which p is nonpositive, for any p G (1, 2) U (2, 2*). If d > 3 and p = 2*, it follows from ([9]) that = 1 
and /3 is uniquely defined. Consistently, we observe that lim p _»2* (/3+ (p) — /3_(p)) = 0, so that the interval 
shrinks to a point when p approaches the critical value 2*, d > 3. Alternatively, there are values of /3 for 
which a whole interval in terms of p can be covered, as it has been noticed in [36] . in the case of the sphere 
when /3 = 1: in such the method applies for any p e (1, 2) U (2, 2 B ], where 2 J := (2 d 2 + l)/(d - l) 2 . 

Our paper leaves open the issue of deciding under which conditions the interval (A, Ai) is empty or not. 
In other words, are there non constant minimizers of T for some A < Ai? If A* < Ai, rigidity results of 
Theorem [3] and optimality results of Theorem 2] are apparently of different nature, although they coincide 
in the case of the sphere. It is not difficult to produce examples of manifolds for which A* < Ai (take non 
constant curvatures and consider the case p in a neighborhood of 2* when d > 3). However, in this paper 
we are dealing only with sufficient conditions and it is therefore not possible to conclude that criteria for 
rigidity are stronger than criteria for determining the optimal constant in the interpolation inequality ([3} . 
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